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ABSTRACT 

We examine the fraction of early-type (and spiral) galaxies found in groups and clusters 
of galaxies as a function of dark matter halo mass. We use morphological classifications 
from the Galaxy Zoo project matched to halo masses from both the C4 cluster cata- 
logue and the Yang et al (2007) group catalogue. We find that the fraction of early-type 
(or spiral) galaxies remains constant (changing by less than 10%) over three orders of 
magnitude in halo mass (13 ^ log My//i -1 M «C 15.8). This result is insensitive to our 
choice of halo mass measure, from velocity dispersions or summed optical luminosity. 
Furthermore, we consider the morphology-halo mass relations in bins of galaxy stellar 
mass M* , and find that while the trend of constant fraction remains unchanged, the 
early-type fraction amongst the most massive galaxies (11^ log M*//i M© 12) is 
a factor of three greater than lower mass galaxies (10^ log M*//i _1 M Q ^10.7). We 
compare our observational results with those of simulations presented in De Lucia et 
al (2011), as well as previous observational analyses, and semi-analytic bulge (or disc) 
dominated galaxies from the Millennium Simulation. We find the simulations recover 
similar trends as observed, but may over-predict the abundances of the most massive 
bulge dominated (early-type) galaxies. Our results suggest that most morphological 
transformation is happening on the group scale before groups merge into massive clus- 
ters. However, we show that within each halo a morphology-density relation remains: 
it is summing the total fraction to a self-similar scaled radius which results in a flat 
morphology-halo mass relationship. 
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1 INTRODUCTION 

One of the central observations that must be explained by 
any model of galaxy evolution is the morphology-density re- 
lation (Dressier 1980): that elliptical (or early-type, bulge- 
dominated) galaxies are more commonly found in the dens- 
est regions of the Universe, whilst spiral (or late-type, disc- 
dominated) galaxies inhabit less dense regions. In hierarchi- 
cal models of galaxy formation, in an expanding universe, 
this observation has been explained by the transformation 
of disc-dominated galaxies into bulge-dominated galaxies via 
heirarchical merging. This basic picture is supported by an 



impressive array of observational evidence, matched to semi- 



et al.|1993||Baugh et al. 


1996; Cole et al.|2000||G6mez et al. 


2003 |Baldry et al.|2004 


Zehavi et al.|2011 


1, although many 


questions remain over the details (e.g., Bower et al. 2010 



has also been a growth of interest in the role internal and/or 
secular evolution can have on the morphological evolution, 
and the growth of spheroidal components in galaxies (e.g., 
Oesch et al.|2010| [Cisternas et al.|2011 l. 



The standard theoretical picture of structure formation 
(gravitational instability of a Gaussian initial density field) 
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predicts that the most massive halos form at the biggest 
peaks of the primordial density field (e.g., |Kaiser|[l984[ ), 
which in turn are more often found in the densest parts of the 
universe. This naively suggests that a morphology-density 
relation will results in a morphology-halo mass correlation. 



In addition, recent observational evidence (e.g., Haas et al 



2012 ) suggests that the most common environmental mea- 



sures used in the literature correlate extremely well with the 
host halo mass (but see Muldrew et al.|2012 who argue that 
nearest neighbour measures are largely independent of halo 
mass), so that the observed "morphology-density" relation 
may in fact be a "morphology-halo mass" relation. Both of 
these arguments seem to suggest that the fraction of bulge- 
dominated galaxies in a given halo should increase with halo 
mass. 



However, De Lucia et al. (20121 recently used a com- 
bination of the Millennium Simulation ( |Springel et al.| 
2005 Boylan-Kolchin et al.|2009a l with two different semi- 



analytic models of galaxy formation, namely DLB07 (De 
Lucia fc Blaizot|[2M7l > and MORGANA ( |Lo Faro etaT 



20091, and found that the fraction of bulge-dominated 



galaxies in both models remains constant, as a function 
of halo mass (Mh), over a narrow cluster mass range of 
14sC logMn/h^Mo s£14.8. 

Observational evidence on the fraction of ellipticals (or 
early-types) in groups and clusters of galaxies as a function 
of halo mass has been limited to relatively small samples due 
to the difficulty in obtaining both reliable halo masses and 
galaxy morphological information. The most recent studies 
suggest contradictory results, probably due to the varying 
halo and galaxy mass scales being probed. Poggianti et al.| 
( [20091 ) agree with the De Lucia et al. (2011) prediction, find- 
ing that the fraction of S0+elliptical (early-type) and spiral 
galaxies remains constant with halo mass (where velocity 
dispersion is used as a proxy for halo mass, limited to the 
range 500 < a v < llOOkms -1 ). This study is based on 
72 low redshift (0.04 < z < 0.07) clusters of galaxies in the 
WINGS survey (Fa sano et al. 12006] ) that possess morpholog- 
ical classifications of galaxies using the photometric package 
MORPHOT, which have been trained on visually classified 
galaxies (see Fasano et al.|2006 for details). 



Calvi et al. (20121 recently studied the morphological 



make-up of galaxy structures (clusters, groups and single 
galaxy halos) as a function of galaxy stellar mass in differ- 
ent halo mass regimes. Their sample contains 176 galaxy 
groups at 0.04 < z < 0.1 as well as isolated, binary and gen- 
eral field galaxies over the same redshift range (all from the 
Padova-Millenium Galaxy Catalogue; Calvi et al.|201l' l. In 
addition, they add galaxies in 21 clusters from the WINGS 
survey ( jFasano et al.|2006| ). This sample spans a halo mass 
range of 12 5j log Mh /h~ Mq ^15, and is limited in galaxy 



Poggianti 



stellar mass (M*) of M* > IO^ ^/i^Mq. As in_ 
|et al.| ( |2009| , the morphologies come from the MORPHOT 
package. They find a smooth increase/ decline in the frac- 
tions of Es-SOs/late type galaxies going from single galaxies, 
to binaries, to groups (i.e. as halo mass increases), which 
does not contradict the predictions of De Lucia et al. (2011), 
since the latter only examine the group/cluster mass scale. 

We extend the above analyses by using an order of mag- 
nitude more groups and clusters of galaxies, and by prob- 
ing a range of (group and cluster) halo masses covering 
13 ^ log Mh //i _1 Mq ^15.8. This paper performs a sin- 



gle, consistent study, over a range of halo masses including 
both groups and clusters, and additionaly considers different 
stellar mass ranges, allowing us to address concerns about 
comparing studies with different halo and galaxy mass se- 
lections. 

To calculate distances, we assume a flat ACDM cos- 
mology with n m ,£l A ,Ho = (0.3, 0.7, 70kms~ 1 Mpc -1 ) and 
h = Ho/100. 

Throughout this paper, we loosely refer to both 
groups and clusters of galaxies as just "clusters" . We also 
use the terms "spiral" and "early-type" galaxy to refer 
to the Galaxy Zoo (GZ) "clean" classifications (as dis- 
cussed below), while the terms "bulge-dominated" or "disc- 
dominated" refer to galaxies in the simulations divided by 
their bulge-disc ratio. We will associate the observed early- 
type class with bulge-dominated galaxies in the simula- 
tions, and likewise the GZ spiral classifications with disc- 
dominated galaxies in the simulations. We remind the reader 
that such an association may have problems as these corre- 
lations are not exact, but should be sufficient to understand 
the broad relationship between observations and theory. 



2 OBSERVATIONAL AND SIMULATED DATA 

All of the galaxies used in this study were drawn from Sloan 



Digital Sky Survey (see York et al.||2000 Gunn et al71|2006 



Smith et al. 2002| and references therein) Data Release 7 
( Abazajian et al.||2009| hereafter SDSS DR7). We made a 
volume limited sample from this catalogue using cuts on 
the k-corrected absolute r band magnitudes M r , taken from 
the PhotoZ table in CasJob^ lOyaizu et al. 2008) with 
-30 < M r < -20 and redshift of z < 0.08. 

We use Galaxy Zoo (GZ) classification probabilities to 
identify the galaxies as either early-type or spiral 



classified). The GZ project ( |Lintott et al.|2008 



201 



used 

an Internet tool to enable citizen scientists to visually clas- 
sify galaxies. The original version of GZ presented users 
with galaxies from the SDSS Main Galaxy Sample (using 
all galaxies which made the spectroscopic targeting criteria 
for the DR7 release) and asked them to be classified as ei- 
ther spiral or early- type (smooth). Each galaxy in this sam- 
ple was classified by a median of 39 citizen scientists (with a 
minimum of 20). The raw results were de-biased (e.g., for the 
effect of higher redshift galaxies appearing smoother as the 
morphological structure becomes blurred), and compared to 
a subset of expert classifiers in Bamford et al.|(|2009). The 



data are now publicly available ( Lintott et al.|20 11p 



When selecting the early-type or spiral galaxy samples, 
we use the de-biased morphological results from |Bamford] 
( 2009 1 who used the citizen scientist classifications to 



et al. 



assign each galaxy with a probability of being an early-type 
galaxy (elliptical + SO) p c i, or a spiral galaxy, p sp . They 
showed that few SO's (predominantly edge-on galaxies) were 
classified as spiral galaxies (i.e., with high p ap ); most SOs 
identified by the experts have high p e i within GZ. We follow 
many of the G Z studies (e.g.,|Bamford et al. 2009 , Schawin- 



ski et al-pOlO] |Masters et al.|2010|?| |2011| ) in using the GZ 



1 http://casjobs.sdss.org 

2 http://www.galaxyzoo.org 

3 |http : //data . galaxyzoo ■ org| 
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"clean" sample, which uses a cut of p sp ,p i > 0.8 to select 
GZ galaxies and define the early- type/spiral separation used 
herein. It is important to remember in interpreting our result 
that this selection results in large numbers of "unclassified" 
galaxies (with both p sp ,p c i < 0.8) which will include the 
irregular galaxies which are included in blue/late- type sam- 
ples of other studies. However, we note that our results are 
largely insensitive to the probability threshold (and there- 
fore the number of unclassified galaxies) , which we demon- 
strate by using the full probabilities in a weighted count (see 
Section 3.4). 

To enable a direct comparison with the |De Lucia et ah] 



(20121 simulations, we additionally impose a stellar mass 
M* , cut of M* > lO lo ft -1 M . The stellar masses were esti- 
mated from the shape of the optical (SDSS fibre) spectrum 



using the publicly available VESPA code (Tojeiro et al 



2007 20091. At the redshifts of the GZ galaxies the SDSS 



fibre captures only a small fraction of the light of the galaxy. 
The VESPA code applies aperture corrections which should 
result in unbiased total stellar masses on average (since as 
shown by Glaz ebrook et al.|2003[ ) for a large sample of SDSS 
galaxies the average fibre colours are very close to the aver- 
age total colours), although the error on an individual galaxy 
stellar mass may be quite large, and this may be a partic- 
ular issue for face-on spiral galaxies which have significant 
colour gradients. To examine this effect we additionally used 



the stellar mass estimates of Baldry et al.| |2004 ) , based on 
integrated SDSS magnitudes and colours. We find that our 
results (in Section 3.2) are insensitive to the choice of stellar 
mass estimate, so conclude that aperture bias errors do not 
critically affect our results. 

Our final galaxy sample after all these cuts consists of 
N = 85364 galaxies, of which TV = 28602 (33.5%) have 
p sp > 0.8 and N = 11571 (13.5%) have p el > 0.8. 

To obtain halo masses, we used the publicly available 
C4 cluster catalogue ( Miller et al.| j'2005) which was applied 
to the SDSS DR5 spectroscopic galaxy sample ( Adclman- 
McCarthy et al.||2007| ). The C4 algorithm identified galaxy 
clusters as over-densities in seven dimensional position (3-D) 
and colour (4-D) space. Unlike many other cluster finding 
algorithms, the C4 cluster catalogue is not biased against 
finding clusters or groups of blue galaxies, and is over 96% 
pure for clusters with an optical (r-band) luminosity of 
L r > 2x 1O 11 L0. However, the sample is only 100% complete 
above halo masses of Mr = 5 x W 14 h~ 1 MQ, and drops dra- 
matically to < 30% complete below M H = 5 x 1O 13 /i _1 M . 

We use only those C4 clusters with a measured virial 
radius, and a halo mass of Mh > lO 13 ft _1 M0. We con- 
verted the virial radius and redshift to an angle subtended 
by the cluster on the sky #vm, and the velocity dispersion 
to a redshift error A(jz). We allocated cluster membership 
to those galaxies which were within ±2A(z), and separated 
by an angle less than #vir from the cluster centre. In total 
N — 10101 galaxies reside within N = 505 C4 cluster halos, 
of which 1848 (18%) are early-type galaxies (p sp > 0.8) and 
2161 (21%) are spiral galaxies (p e i > 0.8). 

We also identify galaxies in halos with masses estimated 
by Y07. These authors estimate the host halo masses of the 
SDSS DR4 galaxy sample ( |Adelman-McCarthy et al.||2006| | 
by iteratively determining the group membership based on 
a luminosity-scaled radius, and then estimate the halo mass 
from a characteristic luminosity, or a characteristic stellar 



mass, of the group or cluster. The mass-to-luminosity re- 
lation was calibrated from simulations ( |Yang et al.| |2004). 
In total, TV = 16567 galaxies in our GZ sample were found 
to reside in N = 2140 Y07 halos (with a mass of Mh > 
lO 13 ft _1 M ), of which 3355 (20%) are early-type galaxies 
(pel > 0.8) and 3908 (24%) are spiral galaxies (p sp > 0.8). 

In this study, we perform two separate analyses, first 
looking at the early-type (spiral) fraction among all galaxies 
in our GZ sample associated with a C4 cluster, and sepa- 
rately then performing the same analysis using all galaxies 
in the sample having a host halo mass estimated by Y07. 

We comment here on the possible differences between 
the C4 and Y07 halo catalogues. For those galaxies with 
both a measured C4 cluster mass, and a Y07 group mass, we 
find a systematic offset of log M H /h~ 1 M Q = 0.5±1.0. This is 
due to the differences in the C4 and Y07 cluster-finding algo- 
rithms resulting in single C4 clusters being "shredded" into 
multiple, smaller, Y07 groups. It is therefore difficult to per- 
form a meaningful one-to-one comparison between C4 and 
Y07 halo masses as C4 halos can be associated with up to 
eleven Y07 groups, and on average, there are (N) = 2.9±3.5, 
Y07 groups (within the C4 virial radius). In summary, it ap- 
pears C4 characterises the largest halo masses, while Y07 is 
optimal for characterising groups of galaxies and substruc- 
tures within larger (C4) systems. A detailed comparison of 
the two halo measurements is beyond the scope of this paper, 
but it is instructive to ensure our results are not dependent 
on the halo measurements we use. 

To compare our data with models of galaxy and struc- 
ture formation, we extracted 960, 000 model galaxies from 
the Millennium Simulation ( |Springel et al.||2005| |Boylan-| 
Kolchin et al.|2009b| , specifically selecting galaxies from the 
DeLucia2006A-SD SS2MASS and DeLucia2006a ( |De Lu-| 
cia & Blaizot 20071 tables. We chose simulated galaxies to 
most resemble those of the GZ galaxies by placing the follow- 
ing constraints on absolute r-band magnitude of Mr ^ —20, 
redshift of z sC 0.08, host halo mass of Mh > 10 13 /i~ 1 M Q , 
and on the galaxy stellar mass of M* ^ lO lo ft _1 M0. We 
identify likely early-type galaxies as those with a bulge stel- 
lar mass to total stellar mass ratio of B/T 0.9 resulting in 
270,000 galaxies (hereafter "bulge-dominated"). Likewise, 
likely disc-dominated (late-type or spiral) galaxies are those 
galaxies with a bulge stellar mass to total stellar mass ratio 
of B/T sC 0.4 (250, 000). Our sample contains 44,000 clusters 
of mass greater than Mh sS 1O 13 /i _1 M0, and 8681 clusters 
within the mass range probed by De Lucia et al. 2011. Our 
sample is much larger because those authors restricted their 
sample to the same number of clusters (100) as available 
from the MORGANA models, which used a smaller simula- 
tion than the Millenium Simulation (see Lo Faro et al.|2009 
for details). 



3 RESULTS 

The fraction of elliptical (or more precisely bulge- 
dominated) galaxies in simulations has recently been shown 
to be constant over the ran ge of cluster halo ma sses of 



14< logMH//i _1 M s£14.8 (De Lucia et al. 20121. This 



result was obtained by identifying bulge-dominated model 
galaxies from an N-body simulation. De Lucia et al. (2011) 
selected 100 dark matter cluster halos and followed the 
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Figure 1. The fractions of early-type (thick red lines) and spi- 
ral galaxies (thin blue lines) per halo with respect to the total 
number of galaxies per halo, as a function of halo mass. We show 
the Y07 halo mass estimates as the solid lines, and the C4 cluster 
halo estimates using the dashed lines. The error bars are Pois- 
son errors. The grey shaded region (labelled "DL Sims") indi- 
cates a prediction for the early-type galaxy fraction obtained by 
simulations (Dc Lucia et al. 2011) which span the mass range 
14< logM H /h- 1 M e <14.8. 



merger histories of their component galactic size dark matter 
sub-halos. They used sophisticated semi-analytic modelling 
of the bulge and disc, and gas cooling flows, to determine the 
present day (z = 0) bulge to disc ratio of the galaxies, and 
identified a galaxy as an early-type if its bulge mass/total 
mass B/T > 0.9 (also our selection above). These results re- 
produced those of previous observations by Poggianti ct al. 
(2009), which were based on the automated morphological 



classifications of a few hundred galaxies in 77 galaxy groups 
and clusters. 



3.1 The dependence of morphology on halo mass 

Our main result is that we agree that the fraction of early- 
type galaxies in a halo is insensitive to halo mass. We show 
this is in Fig. [T] where we plot the fraction of early-type 
(and spiral) galaxies in a halo (as a function of halo mass 
for both the C4 and Y07 halo catalogues). The error bars are 
Poisson counting errors (a good approximation of the errors 
on fractional quantities in large samples which are not close 
to either zero or 100%, iCameron et~al]|2010||. We als o plot 
the result of the simulations by De Lucia et al. (2012 I using 



a constant value for the elliptical fraction, Je = 0.24 ±0.07, 
which is the average result of the two models presented by 
them, and our estimate of the la dispersion in their models 
(from their Fig. 1). The shaded region indicates both the 
range of halo mass in their model and our estimate of the 
dispersion. 

We remind the reader that the fraction of early-type and 
spiral galaxies from GZ, add up to less than one; some galax- 
ies have neither p c i > 0.8 or p sp > 0.8. We find that the frac- 



tion of early-type galaxies is independent of halo mass for 
the C4 catalogue halo mass estimates. This result has been 
seen before in Fig. 13 of| Bamford et al.| ( |2009 [) using the same 



GZ classifications and the C4 cluster catalogue. We extend 
this previous analysis by also showing the independence on 
halo mass also holds using the Y07 halo masses. In addition, 
the fraction of spiral galaxies is shown to decline smoothly 
from 0.3 to 0.2 over the same mass range (again using either 
Y07 or C4 clusters catalogues), and that this process begins 
from halos with masses as small as 131og Mh /h~ 1 MQ. 

3.2 The dependence of morphology on stellar 
mass 

We next examine the fraction of early-type and spiral galax- 
ies as a function of both halo mass and galaxy stellar mass. 
We divide the galaxy sample into bins of total stellar mass 
M* and we compute the fraction of early-type, and spiral, 
galaxies in each bin with respect to the total number of 
galaxies, as a function of halo mass (Mh)- In the top panel 
of Fig. [2] we show these results using the Y07 halo catalogue, 
and in the bottom panel we show the results using the C4 
cluster catalogue. The width of stellar mass bins were chosen 
to contain approximately equal numbers of galaxies. 

In Fig. [2j we see that the fraction of early-type galaxies 
in a halo remains constant as a function of halo mass (within 
the errors) within each stellar mass bin. However, we see that 
the fraction of early-type galaxies per halo increasing by a 
factor of 2 to 3 as we increase the galaxy stellar mass over 
the range 10. 5< logM* jh~ Mq <11.5. It is however unclear 
how much this affect is due to more massive, and therefore 
luminous, galaxies being easier to classify, so the fraction of 
unclassified GZ galaxies which should be classified as early- 
type may drop in the high stellar mass bins. We note that 
the fraction of unclassified galaxies drops as a function of 
increasing stellar mass. 

In contrast, we see that the fraction of spiral galaxies 
in a halo is largely unaffected by either stellar mass, or halo 
mass. We note that we have also examined the effect of split- 
ting the galaxy sample into bins of total galaxy metallicity 
(as derived by VESPA), but find no correlation between 
metallicity and changes to the fraction of either early-type, 
or spiral galaxies. 



Recently Calvi et aLl (|2012|) also found that the fraction 



of elliptical (and the combination elliptical+SO) galaxies in- 
creases as a function of increasing stellar mass, but they find 
that the fraction of late-type galaxies decreases as a function 
of stellar mass. We note that the late-type galaxies of |Calvi| 
et al. ( 2012 I are defined as containing both spiral and irregu- 



lar morphological types, whereas our work has concentrated 
on specifically spiral galaxies from the GZ classifications. 
Particularly at low stellar masses, most late-types may be 
irregulars (e.g., [Binggeli et al.||1988| which would likely be 
"unclassified" in GZ. 

3.3 Comparisons with simulations 

We next present the fraction of bulge-dominated (or early- 
type) galaxies, and disc-dominated (or spiral) galaxies, 
drawn from the Millennium Simulations, as a function of 
halo masses. In the upper panel of Fig. [3] we show the ra- 
tio of the number of bulge (disc) dominated galaxies at each 
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Figure 2. The morphological fractions of galaxy as a function 
of halo mass using the Y07 (C4) catalogue, separated into sub- 
samples of stellar mass in the upper (lower) panel. Here the early- 
type galaxies are shown by the red solid lines, and the spiral 
galaxies by the blue dot-dashed lines. The increasing line thick- 
ness corresponds to increasing bins of stellar mass, and errors are 
Poisson, of which we show every third (all) for clarity. 



halo mass, to the total number of galaxies at each halo mass. 
We again plot the results of the more detailed semi-analytic 



simulations from De Lucia et al. (2012 1 as the shaded region. 
In the lower panel, we split the galaxy samples into bins of 
total stellar mass, with the thicker lines depicting bins of 
increasing stellar mass. We show the error bars as twice the 
Poisson error so they are visible. 

Our results show that the fraction of bulge-dominated 
(or early-type) and disc-dominated (or spiral) galaxies in the 
Millennium Simulation is relatively constant as a function 
of halo mass above logM_tf//i _1 M0 = 13.2, changing by 
at most 10% over two orders of magnitude in Mb (Fig. 
3). We find that the bulge-dominated galaxy fraction is in 
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Figure 3. The fraction of bulge-dominated (early-type) galax- 
ies (red solid lines), and disc-dominated galaxies (blue dot-dashed 
lines) as a function of halo mass, drawn from the Millennium sim- 
ulations. The upper panel shows the ratio of the number of bulge 
(disc) dominated galaxies at each halo mass, to the total number 
of galaxies at each halo mass. We plot the results of the more 
detailed semi-analytic simulations from De Lucia et al. (2011) as 
a shaded region. In the lower panel, we divide the galaxy samples 
into bins of total stellar mass, the thicker lines correspond to bins 
of greater stellar mass as indicated in the key. All error bars are 
twice Poisson to help visualise their size. 



good agreement with the results presented in |De Lucia et al7| 
( |2012[ ), who used a smaller sample size. With these results 
however, we can extend the De Lucia et al. (2011) findings 
over a wider range of halo masses now showing that the 
fraction of elliptical (or more correctly, bulge-dominated) 
simulated galaxies are independent of host halo mass over 
three orders of magnitude. 

The lower panel of Fig. [3] shows how the fractions of 
bulge-dominated (early-type) model galaxies change as a 
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function of halo mass for galaxies in different stellar mass 
bins. We find that the fraction of the most massive bulge- 
dominated galaxies increases by a factor of 3-6 compared 
with the lower mass bulge-dominated galaxies. These trends 
are in qualitative agreement with the data (as seen in both 
panels of Fig.[2|. We do however find that the fraction of the 
most massive bulge-dominated (early-type) galaxies is 50% 
higher than found by observations. The simulation seem to 
over produce the number of massive bulge-dominated (early- 
type) galaxies. However, as noted above, the large fraction of 
unclassified galaxies in GZ (i.e., neither early-type or spiral) 
may introduce an observational bias, because more massive 
galaxies may be easier to classify. We note here that the to- 
tal fraction of bulge-dominated simulated galaxies (28%) in 
massive halos is also greater than the total number of GZ 
early- type galaxies in massive halos (18% for C4, 20% for 
Y07). 

The disc-dominated model galaxies in the top panel 
of Fig. [3] show a fractional change of less than 10% over 
13^ logM^f//i _1 MQ ^15 in agreement with our observa- 
tions. The same mild trend is still apparent if we divide 
the disc-dominated galaxies by stellar mass (see the lower 
panel of Fig. [3|. However, contrary to the data, the total 
abundances of disc-dominated galaxies changes as the stel- 
lar mass increases. This could be real or could be caused by 
the loose comparison of disc-dominated simulated galaxies 
to observed spiral galaxies, or the relatively large errors in 
stellar mass estimates. 



3.4 Using the full GZ probabilities 

Rather than imposing a strict morphological probability cut 
on our galaxy sample (which results in a significant frac- 
tion of unclassified GZ galaxies), in this section we instead 
use the full galaxy sample, weighting each galaxy by the 
morphological probability assigned by GZ. Thus, we calcu- 
late the sum of the probabilities divided by the total num- 
ber of galaxies per halo, as a function of halo mass, i.e., 
fa = E Pel (Mh)MMh) and / sp = E p sp (M H )/N(M H ). 
In Fig. [1] we show the outcome of this analysis, and use the 
same line styles and colours as in Fig. [TJ In agreement with 
Fig.[TJ we again observe that the fraction of early- type (and 
spiral) galaxies is independent of halo mass. 

Considering in more detail the comparison between Fig. 
[4] to Fig. 1 we see that the fraction of both early-type and 
spiral galaxies has increased: as expected since we are now 
including all galaxies, not only a subset with p c i,p sp > 0.8. 
We see that the morphological fractions for both the C4 halo 
masses and Y07 halo masses are consistent with each other, 
differing by only 0.35 sigma for the early-type galaxies (and 

0. 3 sigma for the spiral galaxies). The fractions of early-type 
and spiral galaxies remain relatively flat as a function of the 
halo mass, when using the C4 catalogue, and change by less 
than 20% for the Y07 halo masses. Both early-type galaxy 
distributions are consistent with a constant fraction of 0.46, 

1. e., are independent of halo mass, with a \ 2 value of 0.5 
(0.6) for the Y07 (C4) catalogues. 

We therefore conclude that our results are relatively in- 
sensitive to the choice of classification probability threshold 
used to identify early-type and spiral galaxies from GZ. 




13,0 13.5 14.0 14,5 15,0 15.5 
Halo Mass [log 1Q M e ] 

Figure 4. The probability weighted fractions (see text) of early- 
type (thick red lines) and spiral galaxies (thin blue lines) with 
respect to the total number of galaxies per halo, as a function 
of halo mass. We show the Y07 halo mass estimates by the solid 
lines, and the C4 cluster halo estimates by the dashed lines. Pois- 
son error bars are shown. 



4 DISCUSSION 

The most intriguing result of De Lucia et al. (2011) and our 
above result, is that the fraction of early-type galaxies is 
constant with halo mass, if we integrate the total number 
of galaxies in the cluster out to the virial radius. This may 
appear to be in contradiction with the strong morphology- 



density relation (e.g., seen in Dressier (19801, and using GZ 
by Bamford et al. 2009 and |Skibba et al.| ( |2009[ )), which finds 
that there are more early-type galaxies in denser regions of 
the Universe. 

We can explicitly examine the morphology-density re- 
lation within the cluster virial radius using C4 catalogue. 
We examine the fractions of early-type galaxies internal to 
a halo as a function of the radial distance from the centre 
of each cluster (which we scale to the virial radius). We cal- 
culate the fraction of early-type galaxies as a function of 
angular separation from the cluster centre 9, normalised to 
the angular extent of the cluster on the sky 8vir- To increase 
signal-to-noise we stack clusters in bins of halo mass. In Fig. 
[5]we show the fraction of early-type galaxies in each annulus 
divided by the area of the scaled annulus, as a function of 
scaled radius for different bins of halo mass. 

We find that this scaled fraction of early-type galax- 
ies increases towards the centre of the cluster, reproducing 
the strong morphology-density relation seen previously. The 
relationship is insensitive to the halo mass bin, due to the 
scaling by virial radius, effectively removing the mass depen- 
dence through the mass radius relation. We have examined 
the sensitivity of this result by using an alternative mass- 



radius relationship (from Johnston et al. 20071, and find 
the observed morphology-density relationship to remain un- 
changed. 

So, when we consider the total fraction of early-types in 



© 2010 RAS, MNRAS 000, [TV?? 



Elliptical galaxy fractions in groups and clusters. 



7 



1.4 - 

1.2 
1.0 - 

0.8 

0.6 



° 0.4 - 

0.2 

0.G 



13.0<log M H /h"'M s <14.0 
14.0<log M H /h"'M s <15.0 
15.0<log M H /h"'M s <1S.O 



I 

I \ 
\ i 
t 

I 

\ \ 

\ \ 



0.0 0,2 0.4 0.6 

Scaled angle 8/9 m 



1 ,0 



Figure 5. The fraction of early-type galaxy as a function of 
scaled angular separation from the centre of C4 halos. Fractions 
are divided by the volume contained in the scaled annulus to 
allow comparison between different sized halos. The increasing 
line thickness correspond to increasing bins of halo mass, and 
the errors (are Poisson) are shown for a representative sample of 
points. 



a halo out to the virial radius we find no dependence on halo 
mass; but looking inside the halo we recover the morphology- 
density relation. Put together this demonstrates the self- 
similar nature of the number density profile of galaxies in 
massive halos; by summing the fraction out to the virial ra- 
dius we always end up with the same total fraction of early- 
types, independent of the halo mass. This is not in contrast 
to any morphology-density relation because by definition the 
average density out to the virial radius of a group or clus- 
ter (or any virialised structure in the Universe at a given 
redshift) is the same independent of the total mass of the 
structure (as explained in detail in Poggianti et al. 2010). 
Our results support and extend the findings of Poggianti 
et al. (20101, who perform a joint analysis of simulations 



drawn from the Millennium simulation, and observations of 



227 Abell clusters with SDSS spectra (sec Poggianti et al 



2006 for details), and find a constant number density of 



star-forming galaxies as a function of halo mass. 



5 CONCLUSIONS 

Using a sample of galaxies with visual morphological clas- 
sifications from the Galaxy Zoo project (Lintott et al. 
2008,2011), combined with halo mass estimates from the 
Y07 group catalogue ( Yang et al.|2007 l and C4 cluster cat- 
alogue (Miller et al. 2005) we find: 



• The fraction of early-type galaxies in a given halo 
(defined as having p e i > 0.8 from GZ classifications) re- 
mains constant over a large range of halo mass (13^ 
logM_tf//i _1 M0 ^15.8). This observation holds using either 
halo masses from the Y07 group catalogue or the C4 cluster 
catalogue (also see Figure 13 of |Bamford et al.|2009[ ), and is 



insensitive to the exact details of the classification probabil- 
ity used to identify early-types in GZ. As shown in [j4j this is 
not in contradiction with the morphology-density relation. 

• The fraction of spiral galaxies in a halo (defined as p sp > 
0.8 from GZ classifications) evolves slowly (< 10% for both 
the C4 and Y07 clusters) as a function of halo mass. 

We next divided the galaxy sample into bins of galaxy 
stellar mass M* , and find constant fractions of early- 
type galaxies (or slowly changing fractions of spirals) as 
a function of halo mass in each stellar mass bin. How- 
ever, we find a larger fraction of early-type galaxies (by a 
factor of three) amongst the most massive galaxies (11^ 
log M *//i -1 M0 ^12) in a halo than are seen amongst lower 
mass galaxies (10sC log M7ft _1 M ^10.7). The spiral frac- 
tion in group and clusters size halos appears to be insensitive 
to the stellar mass range of the galaxies considered. 

We compare the above observational results to the |De| 
|Lucia et al.| ( |2012| simulations as well as galaxies and ha- 
los selected from the Millennium Simulation, using the same 
redshift, mass, and absolute magnitude selection criteria as 
the data. We identify the early-type galaxies in the MS as 
those galaxies which are bulge-dominated, with a bulge stel- 
lar mass to total stellar mass ratio of B/T > 0.9. 

Our main results from this analysis are: 

• The fraction of bulge-dominated (early-type) galaxies 
in the simulations is constant (to within 5%) as a function 
of halo mass over the range 13^ log Mh /h~ 1 MQ ^15.6. This 
result is in good agreement with the observations above, and 



with the results of De Lucia et al. ( 2012 1, who used a smaller 



sample of clusters to select bulge-dominated (early-type 
galaxies) in cluster halos (over a smaller mass range of 14?C 
log Afj^//i _1 M© 5C14.8). We can extend their prediction over 
a wider range of halo mass (13^ log Mh /1i~ 1 Mq ^15.6) . 

• The fraction of disc-dominated galaxies in the simula- 
tions drops from 0.3 to 0.2 as the halo mass increases from 
13^ logM/f//i -1 M0 ^15.6, and that this process begins 
from halos with masses as small as 131ogM_f///i _1 M0. These 
trends are observed for the GZ spiral galaxies in both the 
C4 and Y07 halos. 

We again divided the model galaxy sample into bins of 
galaxy stellar mass M* and we found that bulge-dominated 
galaxies in the simulations are more common amongst 
the most massive galaxies (similar to our observed result) 
but the difference between the fraction of bulge-dominated 
(early-type) galaxies amongst the most massive and least 
massive galaxies is ~ 50% greater than we observed. This 
result may, however, be explained by the large fraction of 
unclassified galaxies in the GZ sample, as more luminous, 
massive galaxies in our volume limited sample, are larger, 
and easier, to classify morphologically. The disc-dominated 
galaxies in the simulations are more common (by a factor of 
3) amongst the least massive galaxies in the halos than they 
are amongst the most massive. This is contrary to what we 
observed (that in group and cluster sized halos there was no 
change in GZ spiral fraction with stellar mass), but agrees 
with the observations of |Calvi et al.| ( |2012[ > who found that 
the the most massive late-type (disc-dominated) galaxies are 
less abundant, by a factor of two, in clusters, compared with 
groups. 



© 2010 RAS, MNRAS 000, [J?? 



8 Hoyle et al. 



The flatness of the fraction of early-type (or 
SO+elliptical) and spiral galaxies with halo mass has been 
seen before, albeit over a smaller range in mass (or veloc- 
ity dispersion 500 < a v < 1100 kms -1 ) by Poggianti et al. 
p009l) , usin g 72 low redshift (0 04 < z < 0.07) clusters from 
the WINGS dFasano et al. 12006 1 survey. We extend this anal- 



yses by using an order of magnitude more clusters, and by 
probing a far greater range in velocity dispersion (equivalent 
to 100 ^ fj v ^ 4000 kms -1 ). Our results are also in general 



agreement with the recent paper by Calvi et al. (20121, who 
show that the fraction of ellipticals galaxies remains rela- 
tively flat (at ~ 0.3) between groups and clusters. 

Furthermore, after the submission of this paper, 
Wilman fc Erwin|p012| examine the morphology density re- 
lation, using 911 galaxies with morphological classifications 
from RC3 ( de Vaucouleurs et al.|[l991 1 with further checks 
for morphological type using with SDSS imagining, (note 
that 406 of these had a dubious classification after visual 
inspection, of which 165 were later reclassified) within 729 
Y07 halos of mass 11< logM H //i _1 M <14.8. In that study 
most halos have only one galaxy, and they find a strong de- 
pendence of morphology on halo mass, in particular in the 
fraction of SO galaxies, a classification we are unable to probe 
with GZ1. For comparison, our study is larger, and uses 2140 
halos with 3355 (p c i > 0.8) early-type and 3908 (p sp > 0.8) 
spirals galaxies. 

It is fascinating that while there is a strong morphology- 
density relation (seen above and in Dressier et al. 1980, and 
in GZ Bamford et al. 2009 and Skibba et al. 2009) that 
between group and cluster sizes halos we see no (or little) 
change in the early-type (or spiral) galaxy fraction per halo. 
The interpretation of this result must be that the transfor- 
mation between spiral and early-type galaxies is predomi- 
nantly happening on the galaxy group scale, and very mas- 
sive galaxy clusters are consistent (at least in this obser- 
vation) with having been made simply from a hierarchical 
aggregation of galaxy groups with no further need for mor- 
phological transformation of the galaxies inside the groups. 
That processes which transform galaxy morphologies appear 
to happen preferentially on the galaxy group scale has been 
discu ssed before (|Zabludoff fc M ulchacy 1998 . Kodam aet al.| 



2001 Helsdon & Ponman 2003 Kautsch et al. 20081. Our 



observations and extension of the mass range of the semi- 
analytical models of De Lucia et al. (2011) support this basic 
picture, and suggest that it is well incorporated into basic 
semi-analytic models of galaxy formation in an hierarchical 
universe. 
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